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Abstract 

We introduce completely monotonic functions of order r > and 
show that the remainders in asymptotic expansions of the logarithm of 
Barnes double gamma function and Euler's gamma function give rise 
to completely monotonic functions of any positive integer order. 

1 Completely monotonic functions of positive or- 
der 

In this paper it is established that the remainders in asymptotic expansions 
of the logarithm of Barnes double gamma function and Euler's gamma func- 
tion are (up to a sign) completely monotonic functions of order comparable 
with the decay of the remainder in the expansion. 

We first recall some definitions and give some preliminary results. 

A function / : (0, oo) ^ M is called completely monotonic if / has 
derivatives of all orders and satisfies 

> 0, forall x>0 and n = 0,l,2,... (1) 

J. Dubourdieu [3] proved that if a non constant function / is completely 
monotonic then strict inequality holds in ([1]) . See also [6j for a simpler proof 
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of this result. A characterization of completely monotonic functions is given 
by Bernstein's theorem, see |14l p. 161], which states that / is completely 
monotonic if and only if 

/•oo 

f{x)= / e--*dMt), 

where // is a nonnegative measure on [0, oo) such that the integral converges 
for all X > 0. 

Here we are interested in the class of strongly completely monotonic func- 
tions, introduced in [13]. A function / : (0, oo) ^ M is called strongly com- 
pletely monotonic if it has derivatives of all orders and (— 1)" x""*"^ /^"^(x) 
is nonnegative and decreasing on (0, oo) for all n = 0, 1, 2, . . .. (It is clear 
that being strongly completely monotonic is stronger that being completely 
monotonic.) These functions are connected to the important question of 
superadditivity (cf. [13]). 

The following proposition contains a simple characterization of strongly 
completely monotonic functions and its proof is obtained by a direct appli- 
cation of the definitions given above. 

Proposition 1.1 A function f{x) is strongly completely monotonic if and 
only if the function xf{x) is completely monotonic. 

In [13] the authors gave another characterization of strongly completely 
monotonic functions. 

Proposition 1.2 The function f [x) is strongly completely monotonic if and 
only if 

f-OO 

f{x)= / e--'pit)dt, 
Jo 

where p{t) is nonnegative and increasing and the integral converges for all 
x>Q. 

We notice that any right-continuous function p appearing in Proposition [L2] 
can be written as 

p(t) =/x([0,t]), 

where is a Radon measure on [0, oo). 

An extension of Proposition 11.11 and 11.21 is the following. 

Theorem 1.3 Let r he an integer > 1. The function x'' f{x) is completely 
monotonic if and only if 

POO 

f{x)= / e-^'p{t)dt, 
Jo 
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where the integral converges for all x > and where p is r — 1 times differ- 
entiable on [0, oo) with p^^~^\t) = /x([0, t]) for some Radon measure fi and 
p('=)(0) =0 for0<k<r-2. 

Remark 1.4 The conditions on p ensure that p,p', ■ ■ ■ ,p^^~^^ are all non- 
negative. 

Proof. Suppose that f{x) is completely monotonic. Then 
x^f{x)= / e--*d/x(0 = /:(M)(^), 

for some Radon measure n (and where C denotes the Laplace transform). 
Furthermore, it is a fact that 

x-' = -, — / f-^e-^'dt, 

(r - 1)! Jo 



and this yields f{x) = C{p){x) where p is the convolution of these two 
measures on the half line, 



Prom this formula it is easy to check that p has the asserted properties. 
Conversely, if 

p(^-i)(i) = M[0,i]) 

and all derivatives of p at t = up to order r — 2 are zero then we find by 
integration (r — 1 times) 

^^(*) = r^ fit-sy-'di^is). 

[r - 1)! Jo 

This gives 

x'fix) = x'C{p){x) 



= x^L 



r— 1 

S 



(r-1)! 



^dsj {x)L{^ii){x) = C{n){x), 



and therefore x'^ f{x) is completely monotonic. □ 
In the light of these results we formulate the following definition. 
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Definition 1.5 Let r > {}. A function f defined on (0, oo) is said to be 
completely monotonic of order r ifx'''f{x) is completely monotonic. 

According to this definition, completely monotonic functions of order are 
the classical completely monotonic functions, order 1 are the strongly com- 
pletely monotonic functions and so on. 

Theorem 11.31 is the basis of the proofs of Theorems 12. H 13.11 and 14.11 

Remark 1.6 Following |T|, Section 2.9] a fractional integral Ia{ii){t) (for 
a > of a measure on [0, oo) is defined by 

ia{m = T^^ [\t-sr-Uf,is). 

The proof of Theorem \1.3\ actually shows that f is completely monotonic of 
order a > if and only if f is the Laplace transform of a fractional integral 
of a positive Radon measure on [0,oo), that is, 

fix) = —- e--* / {t-sr-U^^{s)dt. 
Jo Jo 



2 Euler's gamma function 

Let us consider the asymptotic expansion of the logarithm of Euler's gamma 
function 

log T{x) = - logx - X + ^ log(27r) + 

where are the Bernoulli numbers. It has been shown in [2] that for all 
n = 0, 1, 2, . . ., the remainder Rn{x) is completely monotonic on (0, oo). We 
strengthen this result in Theorem 12.11 below. 
For i > we write 

7^ = 1 4 + (Itt + (-1)" 
fc=i ^ '' 



4 



where the remainder term Vn{t) is given as (see |12^ p. 64]) 



oo 
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Vn{t) = ^ 



(t2+47r2 k^){2TTk) 



2n ' 



n > 0. 



(4) 



fc=i 



Theorem 2.1 The function Rn{x) defined in ([2]) is completely monotonic 
of order k for n > k for any k > 0. Indeed, 



where r„(t) = t^^Klt) satisfies r^^^O) = OforO<l<k-l and r^n'{t) > 
for t > 0. 

A similar result regarding complete monotonicity of the remainder of an 
asymptotic expansion of a ratio of gamma functions of the form r(x + 
a)/r(x + b) in terms of powers of + c), is given in Frenzen's paper 



where i?2n+i(^) cLfc the Bernoulli polynomials. We observe that this formula 
extends to t = because of the well-known property 



for n > 0. 

In [8] it is proved that Vn is positive, decreasing and satisfies (t^ > 0. 

We extend the last property in the following Lemma |2.3[ 

Lemma 2.3 The function Vn has the following properties: 




Remark 2.2 In Jl^ it is proved that Vn{t) takes also the form 





(i) (t2^V„(t))(') > for alln>0 and all l<j, 

(ii) (t^J-iKW)^'^ > for alln>0 and all l<j-l. 
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Proof. We have 

oo oo ^ 

where 

= 

Here we notice that s^j\x) tends to zero as for x tending to 0. We get, 
by differentiation, 

1 (/).- ... 1 



for large k. We have Sj{x) = ^j(x^), where 



where the series converges uniformly due to the behaviour of s - {t/{2TTk)) 



0(x) = - 



+ X 



and from Lemma 12.41 below it follows that s^j\x) are all positive. This 
proves the first assertion. The second assertion is proved in the same way. 
(Of course the first assertion follows from the second when I < j — 1.) □ 

Lemma 2.4 Let 

rj.n 

in{x) = ——. 

1 + X 

Then £,n\x) > for k < n and x > 0. 

(k) 

Proof. First of all we notice that (0) = for k = 0, . . . ,n — l. To consider 
the n'th derivative, S,nix) is rewritten as follows: 



1 + X 

n 

\n—k 



k=0 ^ ^ 

k=i ^ ^ 



X + 
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so that 



= (-l)"(x + l)-("+i)(-l)"n! = j-^Yy^ > 0. 
Using the relation 

Jo 

recursively for k from n — 1 to 0, we find that the derivatives of ^„ of order 
not exceeding n are positive. □ 
In the next proposition we gather some additional properties of the func- 
tion Vn that are of interest in their own right. 

Proposition 2.5 The function Vn has the following additional properties, 
(i) t^^Vnit) = i-irVoit) + ELi(-l)"~'K-fc(0)i2^-2, n > 0. 
(U) K(0) = (-ir (1^, n>0. 

(ill) K'(0) = -2K+i(0), n>0. 

(iv) Vl[{t) - Vl[{<d) > 0, Vt > 0, n > 0. 

Remark 2.6 The function Vq appearing in this proposition satisfies by defi- 
nition t / {e*^ - 1) = l-t/2+t'^Vo{t), whence Vo{t) = l/t'^{{t/2) coth{t/2)-l). 

Proof. The first assertion is obtained by repeated application of the recursive 
relation 

t^Vn{t) = Vn^l{0)-Vn-l{t), 

see [H]. The remaining assertions are easily obtained by using definition ([3]) 
and relation □ 
Proof of Theorem \2.1\ Using Binet's formula 

log T{x) = (^x - ^\ogx - X + - log(27r) + 
r°° /t t \ e~^* 

and formula ([3]) we see that 

r„(t) = t2«v;(t), n = 0,l,2,... 

It is clear that rn\^) = for k < n — 1. Now use (i) of Lemma [2.31 to prove 
that r^a\t) > for t > 0. The result then follows from Theorem II. 3[ □ 
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Corollary 2.7 For any n > 1 the function 




.n+l 



Y log(27r) 



n 



B2k 



E 



{2k - 1) 2k 



X 



.n-2fc+l 



fc=l 



is completely monotonic on (0,oo). 

Proof. It follows by a combination of Theorem 12.11 with Theorem 11.31 and 



3 Barnes G-function 

We note that a result similar to Theorem 12.11 holds for the remainder in 
an asymptotic expansion (due to C. Ferreira and J. L. Lopez, [H Theorem 
1]) of the logarithm of Barnes G- function. This function is defined as an 
infinite product and satisfies G(l) =1 and G{z + 1) = r{z) G{z). See also [9j 
for details and additional considerations. The remainder in this expansion 
takes the form 



where Vn{t) is as above. 

Theorem 3.1 The remainder (— l)"P„(x) is completely monotonic of order 
k on (0, oo) for n > k + 1. 



Proof Let A„(t) = t^^-i Clearly A}f^(0) = for < n - 1. It follows 

from Lemma 12.31 that ^^(t) > for t > 0, and the result follows from 



4 Barnes double gamma function 

This section is devoted to the investigation of the remainders in an asymp- 
totic expansion due to Ruijsenaars of the logarithm of Barnes double gamma 
function. The expansion is given in terms of generalized Bernoulli polyno- 
mials Bj^ (x), see [U p. 615] or [12l p. 4]. Our investigation is based on 



□ 




Theorem Ol 



□ 
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Ruijesnaars' results and therefore we have found it natural to his termi- 
nology B2^k{x) = Bk {x). We shall futhermore call these polynomials the 
double Bernoulli polynomials. (Ruijsenaars calls B^^kix) multiple Bernoulli 
polynomials.) 

The double Bernoulli polynomials B2^k{x) are defined by 



^ ' fc=0 



and the double Bernoulli numbers B2^k by -B2,fc = ^2,A:(0). 

The asymptotic expansion of the logarithm of Barnes double gamma 
function with both parameters equal to 1, logr2(w) = log r2(ty|l, 1), is 
given as follows: 

\ogT2{w) = log w + ^-B2,0W^ + B2,lW 

+ E - ^)\B2,kW^-^ + i?2,MM, 

A:=3 

where the remainder R2,m has the representation 

Here, Ku- > and M > 2. See dD (3.13) and (3.14)]. 

In [8| and |10] it was shown independently that R2,2n{x) is a 

completely monotonic function. Below it is verified that it is indeed a com- 
pletely monotonic function of order k for n > k + 1. 

We briefly indicate the two different proofs of complete monotonicity. 
By Bernstein's theorem it amounts to showing the positivity of Un{t) for 
t > and any n > 1, where 

since 

(-l)"-^i?2,2n(x) = e--*^ dt. (5) 

In [8] it was shown that 

Unit) = t2"+V„_i(t) + (t2-+iv;(t))' , (6) 
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where Vn{t) is defined in ([3|) and the proof is obtained by showing that 
{t'^^~^^Vn{t)y > 0. In [lOj the proof is based on a contour integration argu- 
ment and the following representation of Pn{t) = Un{t)/t^ in ([5]) is found 

(2n-l) 2t ' 



27rk t2 ^ (27rA;)2 J ' 

This clearly shows the positivity of Un- 

The main result is formulated in the theorem below. 

Theorem 4.1 The remainder {—l)"^^^ R2,2n{x) is completely monotonic of 
order k on (0, oo) for n > k + 1. 

Proof. It follows from ([6]) that 

Pn{t) = ^ = t2"-V„_i(t) + (2n + l)t2"-iy„(t)+t2X(t) 

= r„_i(t) + A„(t) + r;(t), 

where r„(t) = t^" and A„(t) = t^"-! 

Clearly pl^^\o) = for k < n — 2. Since r^\t) > for < A; < n and 
Ai''^(t) > for < A; < n-1, t > 0, it follows from the above that pli'\t) > 
for < k < n — 1 and for all t > and by Theorem 11.31 this completes the 
proof of the Theorem. □ 



Remark 4.2 Theorem \4-l\ can also be obtained directly from the represen- 
tation ([7]). Indeed, since 

.2n-2 Y^/g,,Ni-2n ^tt/c ^ ^ 2 (t/(27rA;))^"~^ 
^ t^ + {27Tky Z^^{2nkYi + (t/{27:k)f' 

^ (t2 + (2vrA;)2)2 A. (2,^^)4 + (t/(2vrfc))^)2 

and 

,2n-2V., (2n - 1) 2t _ ^ 2(2n - l)t (t/(2vrfc))^"-^ 

' 2nk t^ + {2nkY ^ {2nkY l + {t/{2T:k)f 
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we have from ([7]) 



2 °° At 

k=i ^ ' k=i^ ' 

where gn is defined by 

^2n ^2n 



gn{x) = (n + 1/2)— ^ + 



(l + x2)2' 

We put 

hn{x) = n— h 



'l + X ' + 

and have in this way gn{x) = s„(x)/2 + hn{x'^). The positivity of the n'th 

(k) 

derivative of gn clearly follows from the positivity of the derivatives hn for 
k <n. To investigate the derivatives of hn we rewrite it as follows. 



K{x) = — — + 



x + l ' (x + l)2 
n 



^{(-ir + E(:)(^+i)^Mr- 
+(^^|(-ir+n(x+i)(-ir-i + 

+E(3^-+i)'(-ir-'} 

+ ln{x), 



k=2 

-ly 



where In is a polynomial of degree n — 1. Therefore 
Furthermore, hn\o) = for k < n — 1, whence 



/i(f-i)(x) = / hi^\t)dt>0 



for k = 1, . . . , n. This completes a different proof of Theorem \4.1\ 

Acknowledgment. The authors thank Christian Berg for his comments in 
particular regarding Theorem 11.31 
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